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Abstract 

We present a new non perturbative renormalization group for classical simple fluids. The 
theory is built in the Grand Canonical ensemble and in the framework of two equivalent scalar 
field theories as well. The exact mapping between the three renormalization flows is established 
rigorously. In the Grand Canonical ensemble the theory may be seen as an extension of the 
Hierarchical Reference Theory (L. Reatto and A. Parola, Adv. Phys., 44, 211 (1995)) but however 
does not suffer from its shortcomings at subcritical temperatures. In the framework of a new 
canonical field theory of liquid state developed in that aim our construction identifies with the 
effective average action approach developed recently (J. Berges, N. Tetradis, and C. Wetterich, 
Phys. Rep., 363 (2002)). 
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I. INTRODUCTION 



During the last twenty years the Wilson approach^ to the renormalization group (RG) 
has been the subject of a revival in both statistical physics and field theory Two main 
formulations of the non perturbative renormalization group (NPRG) have been developed 
in parallel. In the first one, a continuous realization of the RG transformation of the Hamil- 
tonian Ti^] is made and no expansion is involved with respect to some small parameter 
of the action TCk[<f\- At scale-fc (in momentum space) the high energy modes <p q , q > k, 
have been integrated out in Ti.^ [ip] which requires the implementation of some cut-off of the 
propagator aiming at separating slow (q < k) and fast (q > k) modes. The flow of the action 
is governed by either by the Wilson-Polchinski equation^* 3 - in case of a smooth cut-off or 
the Wegner-Houghton^ equation in case of a sharp cut-off. These equations, due to their 
complexity, call for the use of approximation and/or truncation methods. We refer to the 
recent review of Bagnuls and Bervillier— of a detailed discussion of this version of the NPRG. 
The second formulation, called the effective average action approach, was developed after 
the seminal works of Nicoll and Chang (sharp cut-off version) 6 and Wetterich (smooth cut- 
off version) 7 . This method implements on the effective action -the Gibbs free energy 
of the fast modes $ g , q > k- rather than on Tik the ideas of integration of high-energy 
modes that underlies any RG approach. The flow of results in equations which can be 
solved under the same kind of approximations than those used for the Wilson-Polchinski 
or Wegner-Houghton equations. Recent reviews and lectures on this second approach are 
available 8,9 and will be consulted with profit. It must be stressed that these two approaches 
of the NPRG are of course equivalent as long as the exact equations are considered but may 
yield different results when approximations are introduced at some stag o 10 ' 11 . The effective 
average action method has been applied successfully to compute with accuracy the critical 
exponents of the Ising model non perturbatively for instanc o 8 ' 12 ' 13 and also to tackle with 
other highly non-trivial problems of condensed matter physics such that reaction-diffusion 
processes^, frustrated magnets^ or the Kosterlitz-Thouless transition^, etc. 

In this paper we implement the method of the effective average action for simple clas- 
sical liquids. This can be done either in the framework of a statistical field theory of 
liquid o 17 ! 18 ! 19 ! 20 ! 21 or in the usual Grand Canonical (GC) formalis m 22 ' 23 ' 24 of classical sta- 
tistical mechanics. The first RG approach to the theory of simple fluids and mixtures was 
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proposed more than twenty years ago in a series of outstanding papers by Reatto, Parola and 
co-workers who developed the so-called Hierarchical Reference Theory (HRT) which gives 
an accurate description of the thermodynamics and the structure of a wide class of fluids 
close to and away from the critical region (for a review see reference 2 ^) . It turns out^ that 
HRT is equivalent to the NPRG with sharp cut-off of Nicoll and Chang&. The HRT will be 
obtained here as a highly singular limit of a regular NPRG equation, i.e. as an ultra-sharp 
cut-off limit. We claim that the theory discussed in this paper is free of the flaws which 
affect HRT in the sub-critical region, i.e. that the coexistence and spinodal curves merge 
into a single curv o 26 ' 27 ' 28 . It is likely -and some arguments will be developed in that sense- 
that the singularity of the limit in which HRT is obtained is the cause of its shortcomings 
giving little hope to improve the theory by considering more sophisticated approximations 
and/or truncations such that those proposed recently 2 ^; see however a little known paper of 
Parola^. 

Our article is organized as follows. In section |ll} we present our NPRG theory in the 
framework of the Grand Canonical of classical statistical mechanics. We consider only the 
case of a fluid of hard spheres with additional attractive pair interactions although the 
generalization to multicomponent fluids would be trivial. After some prolegomena aimed at 
resuming the basics notions of the theory of liquids and at fixing our notations (section ((II Aj) ) 
the exact RG flow of the free energy (3Ak and the direct correlation functions are obtained in 
Sections (jHBj) and (jll Cj) respectively. The derivation of the RG flow given here is elementary 
and does not call for the sophisticated mathematical apparatus developed for and applied 
in, deriving the HRT hierarchy. By the way, the genuine HRT flow is obtained as a singular 
limit of our own smooth cut-off NPRG equations in section (|IID|) where we point out that 
some the HRT flow equations for the direct correlation functions are ill-defined. Section (jlllj) 
is devoted to a transposition of the results of Section (JIT} in the framework of the KSSHE 
theory of liquids. This scalar field theory of liquids is obtained after performing a KSSHE 
transformation ^ 1 ^ 2 ^? 4 ^?^ of the GC Boltzmann's factor and was studied intensively these 
past few year a 1 ^! 1 ?^?! 2 ?! 21 . KSSHE theory is a non-canonical field theory in the sense that the 
coupling between the source (the chemical potential for liquids) and the field is non linear 
which makes a direct comparison with Wetterich NPRG difficult. A mapping of KSSHE 
theory (and thus of the GC formalism) with a canonical field theory is however possible and 
discussed at length in Section (jIVj) . Our key result is that the flows of the average affective 
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action Tk of the canonical theory (governed by Wetterich equations) and that of the free 
energy f3Ak (governed by the equations obtained in Section (II)) coincide (cf Section ([IV C)) . 
This mapping makes possible the transposition of all known-properties of the NPRG flows 
for canonical scalar field theories to the theory of liquids. An important consequence of 
this property is that our smooth cut-off NPRG theory for liquids yields a non-vanishing 
compressibility along the coexistence curve and is thus free of the shortcomings of HRT 
along subcritical isotherms (see Section (jIVD)0 . As an illustration of the issues of Sections 
(II, III, IV) we study the case of the Kac mode l 22 i 31 i 37 i 38 i 39 i 40 in Section (V). In this case the 
flow can be solved exactly by two complementary methods. We conclude in Section (VI). 

II. THE RG IN THE GRAND-CANONICAL ENSEMBLE 
A. Prolegomena 

1. The model 

We consider a simple fluid made of identical hard spheres (HS) of diameter a with ad- 
ditional isotropic pair interactions vir^ (r^ = \x{ — Xj\, position of particle "z"). Since 
v (r) is an arbitrary function of r in the core, i.e. for r < cr, one can assume that v{r) has 
been regularized in the core in such a way that its Fourier transform v q is a well behaved 
function of q and that v (0) is a finite quantity. We denote by Q C M. D (D dimension of 
space) the domain occupied by the molecules of the fluid. For convenience Q is chosen to 
be a cube of side L and periodic boundary (PB) conditions are imposed so that the volume 
of Q is V = L D . 

The fluid is at equilibrium in the grand canonical (GC) ensemble, (3 = is the inverse 

temperature (ks Boltzmann's constant), and fi the chemical potential. In addition the 
particles are subject to an external potential i/j(x) and we will denote by v(x) = f3(fi — tp(x)) 
the dimensionless local chemical potential. We stick to notations usually adopted in standard 
textbooks devoted to the theory of liquids (see e.g.—) and thus denote by w{r) = —f3v{r) 
minus the dimensionless pair interaction. Moreover we restrict ourselves to the case of 
attractive interactions, i.e. such that w(q) > for all q. 

In a given GC configuration C = (N;x\ ...xn) the microscopic density of particles at 
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point x reads 



A? 



p{x\C)=Y / S^(x-x i ) , (2-1) 



i=i 



and the grand canonical partition function (GCPF) S [v] which encodes all the physics of 
the model at equilibrium is defined as 2 ^ 

S[u) = Tr[exp(-/W GC [C])] , 

-PHgc = -PVns + \ (p\w\p) + (F|p> , 

oo - 

Tr t---] = E^n IJl - dn- - , (2-2) 



JV=0 



where i = Xi and = d £«. For a given volume V and a given inverse temperature /3, 
H [v] may be considered as a functional of the local chemical potential v{x)^ which we have 
stressed by using a bracket. In equation ()2.2|) /3Vhs [C] denotes the HS contribution to the 
configurational energy and V = v + v$ where i>g — —w(0)/2 is j3 times the self-energy of a 
particle. We have employed here convenient brac-kets notations 

(u\p) = f dl 77(l)p(l) (2.3) 
Jn 

(p\w\p) = /dld2p(l|CMl,2)p(2|C), (2.4) 
Jn 

but, at some places, we will also make use of matricial notations, e.g. 

(pMp)EEp(l).Ml,2)-p(2). 

Finally thermal averages of configurational quantities A [C] in the GC ensemble will be 
noted 

/Ain , Tr [A [C] exp (-gftgg [C])] 

( ^ [CJ)gc " Tr [exp (-pHec ^ ' ( > 



2. Correlation functions 

As a consequence of a long history, the zoology of correlation functions in the theory of 
liquids is quite involved and two families of ordinary and truncated (or connected) density 
correlation functions are considered in the literatur e - 2212 ? 124 . The first family is defined by 
means of functional derivatives of the GCPF with respect to chemical potentials: 

GW[u](l,...,n) = (n^l^) Gc ^ ^ 5[ Mn) ' 
G^ T W)(1 n) = ^ lQ S S M (2 6) 



Our notation emphasizes the fact that the G^[u](l, . . . ,n) are functionals of the local 
chemical potential u(x) and functions of the coordinates (1, . . . , n) = (xi, . . . , x n ). As it is 
well known we hav o 23 i 24 

G^ T [u](l, ...,n) = G<">[i/](1, ,n)-J2U G (m) ' T M(ii, ■••.<»»), (2-7) 

m<n 

where the sum of products is carried out over all possible partitions of the set (1, . . . , n) into 
subsets . . . ,i m ) of cardinal m < n. Of course p[i/](x) = G^ n=l > [u] (x) = G^ n=1 ^ T [z/](x) is 
the local macroscopic density of the fluid. 

The second family of correlation functions is obtained by considering derivatives of the 
GCPF with respect to the activity z(i) = exp(u(i)). These objects will be denoted here by 
the symbol p^ and are defined as 

IN N 

S[v] 5z(l)...5z(n) ' 
,<»>.>](!,..., n) = n^T^ffr- (2-8) 

The ordinary p^ and their connected parts p( n )> T are related by relations similar to (|2.7j) . 
Of course p( n=1 ) = G^ n=1 ^ but, at higher orders, p^ and differ by a host of terms 
involving products of 5 functions and local densities. We need to consider here only the case 
n = 2 for which we have however the simple relation 

G( 2 )(l,2) = p( 2 )(l,2)+p(l)5(l,2)=p(l)p(2)^(l,2)+p(l)5(l,2), 
G< 2 >- T (l, 2) = p(l)p(2)/i(l, 2) + p(l)5(l, 2) , (2.9) 

where g(l,2) (= g(r 12 ) for a homogeneous system) is the usual pair correlation function of 
the theory of liquids, and h(r) — g(r) — 1 tends to zero as r — > oo. 

3. The GC free energy and the direct correlation functions. 

We review here some properties of the GC free energy and the direct correlation functions 
which will be used in the next sections. 
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It can be proved from first principles that the grand potential In S [is] is a convex functional 
of the local chemical potential is(x) at given (3 and V 41 ^ 2 . As a consequence of convexity 
In S [is] is continuous, its functional derivatives exist and are also continuous and it admits 
a Legendre transform (3 A [p] with respect to v defined as 

f3A[p] = sup{(p|z/>-ln~[z/]} . (2.10) 

Although In H [is] can be defined for any well behaved functions is(x) (i.e. with some pedantry 
for is G A/") the supremum in the left hand side of equation (j2.1())) exists only for a restricted 
class 1Z of functions p(x) (for instance if, for some x G Q one has p(x) < then p £ 7Z). 
However if the supremum exists then it is unique and is realized for a unique is*(x) given by 
the solution of the implicit equation: 

5 InE [is] 



(2.11) 



Mi) 

f3 A [p] is then itself a convex functional of the local density p(x) on the convex set 1Z and is 
named the GC or Kohn-Sham free energy. Note that for a homogeneous system f3A [p] and 
the usual Helmholtz free energy f3F(p) of the fluid coincide only in the thermodynamic limit. 
For instance, at temperatures below the critical temperature of the liquid-vapor transition 
and at a finite volume V, (3F(p) will exhibit a van der Waals loop while (3 A [p] will be 
convex (more precisely (3 A [p] is then the convex envelop of f3F(p)) 41 ^ 2 . Since the Legendre 
transform is involutive we have for all v G M 

InSM = sup {{is\p)- (3 A[p]} = (p*\is) ~(3A[p*] , (2.12a) 
pen 

is(l) - ^ 



5p(l) 



(2.12b) 

p=p* 

Let us suppose that In E [is] and (3 A [p] depend in addition on some external parameter, say 
a scalar field a(x). It follows then from the stationary conditions (|2.11j) and f!2.12b|) that^ 



S InE [is] 



Sa(l) 



S[3A[p] 



6a(l) 



(2.13) 



where the partial functional derivatives of InS [is] and (3 A [p] with respect to a(x) must be 
taken at fixed v and p respectively; of course equation (|2.13J) is valid only if v and p are 
related by one of the two stationary conditions ()2.11|) or (|2.12b|) . 
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It was shown in ref.— that in the case of an attractive interaction (i.e. w(q) > for all 
q) we have the rigorous and useful inequality 

(3 A [p] < (3Amf [p] , 
/3A MF [p] = Mrs [p] ~ \ (pHp) - (pWs) ■ (2.14) 

Equation (j2.14j) is valid for any p e 1Z and *4hs [p] is the GC free energy of the HS reference 
system. The subscript "MF" which stands for "mean-field" emphasizes that ^4mf is nothing 
but the van de Waals free energy 2 ^. Note that (3Amf [p] is a priori non convex, notably in 
the two-phase region. 

The n-point direct correlation function is defined in this paper as 

0«W(1 ») = -77^fn. ("S) 

5p(l)...6p(n) 

The definition (J2.15|) of includes an ideal gas contribution which is sometimes discarded 
explicitely by introducing a small cap c^ n ' defined a o 23 ' 24 

( 1 \n n 

C?W(1, . . . , n) = c W(l, . . . , n) + A_J_ ft 6(1, i) . (2. 16) 

It follows from the properties of the Legendre transform that the linear operators G^' T and 
are positive definite and the inverse of each other: 

G (2) ' T (1, 3) . C {2) (3, 2) = -6(1, 2) . (2.17) 

The above relation, when re-expressed in terms of h(l,2) and c(l,2) = c^ 2 ^(l,2) coincides 
with the well known Ornstein-Zernicke (OZ) equation of the theory of liquid o 22 ' 23 ' 24 

h(l, 2) = c(l, 2) + c(l, 3) . p(3) . h(3, 2) . (2.18) 



B. Renormalization group flow for the effective energy 



1. Basic assumptions 

The trace operation which enters the definition ()2.2j) of the GCPF 5 [u] may be seen 
as a sum over all the modes (or collective variables) p q , i.e. the Fourier coefficients of 
the microscopic density p(x) (cf equation (j2.1|) \ These modes are coupled via the HS 
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interactions and the pair potential w(q) as well. Following the ideas which are the corner 
stone of the HRT theory^ we build a one-parameter family of models, called k-systems, 
which are indexed by a scale k in momentum space. The k-system is a HS fluid with pair 
interactions w k (q) which couples the fast modes p q (q > k) as in the original model but not 
the slow modes (g < k). Therefore we ask that w k (q) — w{q) for q > k and w k (q) "C w(q) 
for q < k. Thus, as k decreases, more and more density fluctuations p q are integrated out. 
Denoting provisionally by S k [u] and by f3A k [p] the grand potential and the GC free energy 
of the k-system respectively we require that 

lim (3 A k [p] = f3A[p] , 

k— >0 

lim /3A k \p] = (3Au?[p\ , (2.19) 

k~ +00 

where the MF Kohn-Sham free energy (3Amf [p] has been defined in equation (|2.14|) . Note 
that the limit k — > 00 in equation ()2.19)1 could be problematic and, in practice, an ultra 
violet (UV) cut-off k max ~ 1/cr on vectors k can be introduced if needed. 
In order to implement these ideas let us define 

- • (2 ' 20) 

1 +w(q)R k (q) 

where the regulator R k (q) acts as a smooth IR cut-off for the pair potential and the free 
energy (i.e it suppresses the contributions of the slow modes p q , q < k to (3A k [p] )■ Any 
function of the type Rk{q) = k 2 (l — Q e (q,k)) where Q £ (q,k) is some smooth version of 
the step function Q(q — k) with width w e will obviously do the job. More precisely one 
demands that Q € (q,k) ~ 1 for q > k + e and Q e (q,k) ~ for q < k — e. One could 
think of e (g, k) — (1 + tanh ((q — fc)/e)) /2 for instance. The choice of an analytical form 
for R k {q) should be unimportant in principle but must be considered with some care when 
approximations are introduced at a further stage. Let us summarize the behavior of R k {q) 
in the various regimes of interest: 

Rk(q) > Vg , (2.21a) 
R k {q) w k 2 for q < k , (2.21b) 
R k {q) -> fastly for g -> 00 . (2.21c) 

The positivity condition ()2.21a|) ensures that, at scale-k the pair potential w k (r) is attractive, 
but "less" than w(r), since one has < w k (q) < w(q) Vg. Indeed, as required w k (q) ~ 



w(q) for q > fc and w k (q) ~ for q < fc are consequences of equations (j2.21b|) and (|2.21c|) . 
Finally we note that 

lim.Wk(q) = w(q) , (2.22a) 
k— >o 

lim w k (q) = , (2.22b) 

fc— >oo 

hence, when the flow starts at fc = oo one deals with a fluid of pure hard spheres and, when 
it stops at = one recovers the full interacting system. 

Parola has proposed his own smooth cut-off version of HRT equations 30 , however his 
choice for w k , i.e. w k {q) = w(q) — k 2 ~ v w(q/k) (rj > Fisher's exponent) does not enter the 
class of potentials defined at equation (|2.2(J|) . 

S fc [u] will be defined according to equation (|2.2|) with the replacement w — > w k . Since 
it is the GCPF of a HS fluid with additional pair interactions Wk{r), lnHfc [u] is a convex 
functional of the local potential v(x). Denoting by (3A k [p] its Legendre transform, one has 
the usual couple of relations 

In ~ k [v] = sup { (v\p) - @A' k [p] } W G Af, 
pell L J 

(3A' k [ P ] = sup {(p\u) - lnE k [u]} Vpell. (2.23) 

Recall that f3A k [p] is then a convex functional of the density. However (3A k [p] has not the 
required behavior (|2.19jl for k — > oo and, following Reatto et aL— , one is led to introduce 
the effective average free energy 

(3A k [p] = f3A' k [p] + - ( P \w k -w\p) + (u Si k - u s \p) , (2.24) 

where 1/$, k — —w k (0)/2 is the self energy of particles of the k-system. Now, for k — ► 
one has w k — > w E k — >• H ^> /3^4^ — > f3A' k=0 = (3 A. Moreover, for — >• oo one has 
u>fc — > =^> S fc — > S H s =>• — >■ /^^4hs =^ /^^fe — /?^4mf as required. 

Note that since w k (q) — w(q) < 0, Wq then /3^4.fc [p] is the sum of a convex-up functional 
(i.e. .Afc [p]) and a convex-down quadratic form, hence, in general, it will not be convex 
notably in the two phase region at fc ^ (in particular at fc = oo); convexity will be 
recovered only in the fc — >• limit 8 . 

Before deriving the flow equations for A k [p] let us mention some important properties of 
this flow and precise its initial conditions at fc = oo. Firstly we note that it follows readily 
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from equations (j2.14j) and (|2.24|) that for any scale-k (for fc = 0we already know the result) 
and p G R one has 

(3A k [p\ < (3A UF [p] = PA k=00 [p] , (2.25) 

from which however it would be hazardous (and in fact wrong in general but true in the 
toy model of Section V) to infer that j3A k [p] is a decreasing function of k for any density 
p(x) G TZ, even a uniform one. 

The "true" direct correlation functions at scale-k, noted C k , are of course defined 
according to equation (j2.15|) with the replacement A — > A' k . In addition we also introduce 
effective direct correlation functions as 

^•■■■•^-im^hi- (2 ' 26) 

It follows from equation (|2.24|) that the C k and the C k ^ coincide for n > 3 but that for 
n = 2 one has 

C; {2) [p](l,2) = cf[p](l,2) + w k (l,2) - w(l,2), (2.27) 

from which we infer from the convexity of f3A k [p] that, in k-space 

C ( k 2 \q) + w k (q) - w(q)<0 Vg . (2.28) 

Finally, since (3Ak=oo = ^mf the initial conditions for the are: 

C£L W (1) = 4s } M (1) + u»(l, 2) ■ p(2) + u s , 
C£L\P]{1,2) = C§l[p] (1,2) +w(l,2), 
Ctloo [P] (1, • • • , n) = [p] (1, • • • , n) for n > 3 . (2.29) 

(2) 

Note that the initial condition for Coo [p] (1,2) coincides with the so-called random phase 
approximation (RPA) of the theory of liquid o 1 ^ 2212 ^ . We will assume that the direct corre- 
lation functions of the HS fluid are exactly known which is of course not actually the case; 
however some excellent approximations for the Cgg can be found in the literature 22 . 
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2. Flow equations for In 3 k and (3A k 



It follows readily from the definition ()2.2|) of the GCPF E k [u] (with w — > w k ) that the 
partial derivative d k of In [z/] with respect to scale-k (at fixed v) is given by 



<9 fc In S fc [i/] = 




= ^Jdld2 d k w k (l, 2) p< 2) H (1, 2) , (2.30) 

f2l _____ 

where p k [f](l, 2) is the pair correlation of the k-system as defined in equation (J2.8|) . From 
the stationary property ([2.13)1 with a(x) = k we deduce that 

d k pA' k [p] = ~ J dl d2 d k w k (l, 2) pf [v] (1, 2) , (2.31) 

where p and v are related through the equations (j2.11|) and ()2.12b|) . Injecting now equation 
(12.31)) in the definition (|2.24|) yields the flow equation for the effective free energy, i.e. 

d k (3A k [p] =-ljdl d2 d k w k (l, 2) G[ 2) t [u] (1, 2) , (2.32) 

where we have made use of equation ([2.9)1 . Note that equation ([2.32)1 is somehow reminiscent, 
at least formally, of the charging process of the theory of classical fluids^. Recall that G k 2 ^ T 

i (2\ _____ 

is minus the inverse of C k and, taking into account the relation ([2.27)1 . we get the flow 
equation 

d k f3A k \p] = - J dl d2 d k w k (l,2) {cf [p] + w k - (1,2) . (2.33) 

This beautiful equation which was obtained in a trivial way is exact and thus horribly 
complicated. Mathematically it is a functional partial derivative equation (both f3A k and 

(2) 

C k are functionals of p(x)), needless to say impossible to solve exactly. But two types of 
approximations have been be proposed to solve it, at least numerically; they will be briefly 
evoked in Section VI. 

The one-loop structure of the flow equation ()2.33j) should be more easily comprehended 
if rewritten as 

d k (3A k = -i^Trln(-Cf - w k + _) , (2.34) 

where d k . . . = d k w k d . . . / dw k acts only on the k-dependence of w k and not on C k and 
the trace operation is understood in the usual matricial sense, i.e. as an integral over the 
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position variables 1 and 2. Note that the signs have been chosen in such a way that the 
argument of the In is positive definite. It is the place to mention that it is found that, in 
KSSHE or CV field theories, the one-loop approximation for the "true" GC free energy (3A' k 
of the k-system reads a n 19 ' 20 ' 21 

M (1 - loop) IP) = Mmf, * [p] + ^Tr In (l - w k • G§ T ) , (2.35) 

where @Amf, k [p] is given by equation (j2.14j) with w — > w k . Note that the expression ()2.35|) 
of (3A' k coincides with that obtained in the RPA approximation of the theory of liquids^ 2 -^. 
Taking into account equation ()2.24|) one finds for the average free energy 

M 1-l0 ° P) H = ^mf [p] + ^Tr In (((<z» eg + w) + w k - w) • <7<g - 1 ) , (2.36) 
where we have made use of the OZ equation C^g _1 = — G^l T for the reference HS fluid. 

(2) / (2) ' (2) (2) 

Thus, substituting C^+w (i.e. the RPA approximation for C k since C k mpj^ = C^+w k =>- 
C k (rpa) = ^hs + w ) by the full C k turns the one-loop result into an exact one ! 

For a homogeneous system (3A k [p] = Vfk(p) where (3 times the free energy per unit 
volume fk(p) is a function (not a functional) of p and the flow equation for f k (p) reads as 

2 A Q (?) +Wk{q) - u>(g) 

where J q = J d D q/{2-n) D and Cj?\q) is the Fourier transform of C^(r). Note that the right 
hand side is non singular since the denominator is negative definite (cf equation (J2.28j) \ It 
is a general property that the flow of f k has neither UV or IR singularities, in particular IR 
singularities (near a critical point) are smoothened by k and they build up progressively as 
the scale-k is lowered^. 

One can deduce from equation (J2.33)) an infinite tower of equations for the 
C k [p] (1, . . . , n). This is the subject of next section ()11 C|) . 



C. RG flow of the effective direct correlation functions 



1. Real space 

(n) 

The RG flow for the effective direct correlation function C k is obtained by taking n 
successive functional derivatives of both members of equation ()2.33|) with respect to the 

13 



density p(x). Let us first work out the case n — 1 explicitely. Since (Cjp — w + w k ) 1 = 
—G^' T one has, by taking the functional derivative of equation ()2.33|) with respect to p(l) 

d k c£\l) = \ d k M^l) 5G \l^ ] ■ (2.38) 

It follows from the general result 8 A' 1 (1 ,2) = l') • 5A(l',2') ■ A' 1 (2' ,2) valid for 

any linear operator A(l, 2) and variation S that 

6G %^ ,2) = T (l, 1') • cf (!', 2', 3) ■ G?>' T (2', 2) , (2.39) 

yielding readily 

9^(1) = \ d kWk (l\2).G^ T (l\l).C^\l\2\l).G^ T (2,2) . (2.40) 

The flow equation for C^ n=2 ^ can be obtained by taking the functional derivatives of equation 
()2.40|) with respect to p(2). But we will not do it in that way. In fact, the structure 
of the hierarchy of the flow equations for the C^ n ' is more easy to grasp on the basis of a 
diagrammatic analysis and equation (|2.40J) will be rederived now by using Feynman diagrams 
and simple diagrammatic rules. These diagrams will be built out from two kinds of vertex 



— -O— = d h w k {l,2), (2.41a) 
1 2 

2 V 1 = d k n \l,2,...,n) (n>l) , (2.41b) 

and the propagator 

= Gl 2) ' T (l,2). (2.42) 

1 2 

We can now represent the flow equation ()2.33|) diagrammatically as 

d k (3A k = -~^^} . (2.43) 
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The action of a functional derivative with respect to p(l) on the vertex and the prop- 
agator G k 2 ^' T are given by simple diagrammatic rules: 



1 

(2.44a) 



Sp(n) y y . 1 

' y ' n n n + 1 



3 

(2.44b) 



M3) 1 2 1 2 

where equation (|2.44ajl follows from the very definition of the direct correlation functions 
and equation (j2.44b|) is the diagrammatic representation of equation ()2.39|) . Taking now the 
functional derivative of equation ()2.43|) with respect to p(l) by taking into account the rules 
(glU) one finds that 

d k 1 • = ~ ^^---1 , (2.45) 

which is nothing but the diagrammatic representation of equation ()2.4U|) . Continuating the 
process one finds 



which can be translated algebraically by '2 

d k C (2 \l,2) = ±d k w(a,b)G k 2) ' T (b,c)G k 2) ' T (a,d)cj?\c,d,l,2) 
+ d k w(a,b)G k 2) ' T (a,c)G k 2) ' T (b,d) x 

x ct\c,e,2)G k 2) > T {eJ)Ct\f,d,l). (2.47) 

Flow equations for of higher orders are obtained in the same vein by making use ad 
libitum of the diagrammatic rules (|2.44|) . Some comments are in order. 

• The equation for d k C^ involves inter alias the vertex C^ 1 and C^ k n+2 \ therefore 
the hierarchy never closes. Possible approximations consist in enforcing a closure at 
some order r ^ f ^--~^fe . 
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• Contemplation of equations ()2.45[ I2.46J) and a little thought reveal that the one-loop 
structure is present at each order n and therefore only one integral on internal variables 
survives. 

• All the expressions for the odd d k C^ n+1 ^ include diagrams with at least one odd vertex 
C^ m+1 \ Therefore if at any scale k (in practice k = oo) all the odd c( 2m+1 ) happen 
to vanish they will remain exactly zero at smaller scales k. This property is present 
only for peculiar models of fluids with Z 2 symmetry (for instance some versions of the 
lattice gas with special particle-hole symmetry), but not in general. 



2. Fourier space 

We restrict ourselves to homogeneous systems. Then, as usual, in momentum space, we 
factor out and evaluate the momentum conserving 6 function so n-point correlation functions 
& n \qi, . . . , q n ) are defined only when qi + . . . + q n — 0. More precisely one has for instance 

G {n) ' T (qi, ...,q n ) = <%i + ...q n ) J exp (i {q x xi + ... g n _ix n _i)) x 

x GW> T (x 1 ,...,x n _ 1 ,0), (2.48) 

where 5(q) = (2n) D 5 D (q) and J = f n d D x. In addition, in two-point functions, we solve 
qi = —q2 = q and recognize that they are functions only of q = \q\, e.g. G^' T (q). Feynman 
diagrams in Fourier space will be build from the vertex 

= d k w k (q), (2.49a) 



q q 

'12 x. s qi 

Qn 



_^*C = C k n) (q 1 ,...,q n ) (n>l) , (2.49b) 

and the propagator 



-<? =5(3), T 



G k (<?) = ~ (2) . ; , „ „ ■ (2-50) 



Cr (q) - w(q) + w k (q) 
One has now for the effective average free energy per unit volume (cf equation (|2.43Jl ) 



q 



*A = -Kj 'i/ ffln^l-o ' (2 - 51) 
1 v^y z J iC k (?) - w (q) + w k{q) 



-q 
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One then easily obtains the flows for the low order effective direct correlation functions in 
momentum space. With the convention that the impulsion is conserved at each vertex and 
that all the diagrams drawn below are shorn of their external lines one obtains 



q -q 

~m 10 

d k cW(0) = d k . • = - — 




-q q 



\ j C^{Q,q,-q)G k {q) 2 d k w k {q) , (2.52) 



and 



d k C^{p) = d k I . I 




-q q 




2 



+ 

—p — q~p + q 
C^(p,q, -q, -p)G^ T {q) 2 d k w{q) 
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+ / C (3) (p, q, -p - q)C {3) (-q, -p, p + q)G®> T (q) 2 G {2) > T (p + q)d k w(q)(2.53) 

Jq 



D. HRT as an ultra-sharp cut-off limit 

In HRT the pair potential of the k-system is defined as^ 



w k (q) = w(q) for q > k , 

w k (q) = for q < k , (2.54) 
which, in our theory, corresponds to the ultra-sharp cut-off 

Rk(q) = for q > k , 

R k (q) = oo for q < k . (2.55) 

It is the place to mention a relatively unknown paper of Parola^ where the author proposes 
his own smooth cut-off formulation of HRT theory. Unfortunately his expression for w k 
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does not belong to the class of pair potentials defined at equation ()2.20j) so comparisons are 
hardly possible. Of course in the framework of the genuine HRT some care is needed to 
obtain the flow dkfk from equation (J2.37j) in the singular ultra sharp cut-off limit. In order 
to achieve this aim let us rewrite 

w k (q) = Q e (q,k)w(q) , 
d k w k (q) = -8 e (q,k)w(q) , (2.56) 

where Q e (q,k) and 5 e (q,k) = — d k Q t (q, k) are smoothened versions of the step and Dirac 
distributions respectively. HRT corresponds to the limit e — ► in which Q e (q, k) — > Q(q — k) 
and S t (q, k) — > 5(q — k) . In the case of a homogeneous system the flow of the specific GC 
free energy is thus given by 

4AM = ~ / ™ ^mMl . (2 . 57) 

To extract the limit e — > of equation (|2.57j) we will make use of the "little lemma" of 
Morris 10 which states that, for e — > 

S e (q,k)f(e e (q,k),k)^S(q-k) I dt f(t,q) , (2.58) 

provided that the function f(Q e (q, k), k) is continuous at k = q in the limit e — > 0, which is 
the case here. Applying lemma (|2.58J) to equation (|2.57J) one gets 

l r r 1 , l 



dkfk(p) = -- / w(q)5(q-k) / dt 



q 



l + ^ fc (A;)^J , (2.59) 



V-i_%_ln 

2 (2vr) D 



where So = 2n D ^ 2 /T(D /2) is the surface of the D-dimensional sphere of unit radius and we 
have adopted the notations of Parola and Reatto^ ^Fk{q) = ~^ / C^\q) ■ Equation (|2.59j) is 
nothing but HRT equation for the GC free energy as expected. 

Similarly, for the first equation of the hierarchy in k-space, i.e. equation (J2.52j) 

8 t dP(0) = -lim i / Cf\q, -q, 0)— MM^M _ (2 M) 
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Ci 2 \q)+e e (q,k)w(q)-w(q) 
the application of lemma (|2.58jl yields a well-defined limit 
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which again coincides with the HRT result. 

The case n = 2 is more problematic (and we'll stop here our investigations) since there 
are two contributions I and II to the flow of d k C^ k \p) = 1 + 11 (cf equation ()2.53|) ). The 
e — > limit of the first contribution 



1 



1 = / ct\pi -Pi 9, -q)j—— 

Jq [Ci 2) (q) + Q e (q,k)w(q)-W(q) 

causes no trouble and application of Morris lemma yields the limit 



S t (q,k)w(q) 



[2.62) 



I 



-V- 1 

2 



Cl 4 \p, ~P,k,-k)- 



F h {k) 2 w(k) 



(2.63) 



(2vr) D " ' '1 + T k {k)w{k) 

where the angular measure dVtt denotes the integration over the Euler angles of unit vector 
k = k/k. However the second contribution to d k C k (p) which may be rewritten as 

5 e (q,k)w(q) 



II 



Ck\p, 9, ~P ~ q)C£'(-p, -q,p + q) 



(3), 



Cf\q) + Q e {q^)w{q)-w{q) 



X 



Cf\p + q) + &e(p + q, k)w(p + q) - w(p + q) 



(2.64) 



causes some troubles. Admittedly the first fraction in the RHS of the above equation is OK 
and its limit e — > may be taken safely by applying Morris Lemma but the second one 
shows a discontinuity for || p + q\\ = k. However, for a generic k 7^ 0, the discontinuity 
occurs only in a zero measure subset of the integration domain. Except exactly at k = 
where the integral is not defined, this difficulty would not have any serious consequences as 
claimed by the authors of reference^. Following them we define the (discontinuous) function 
Fk{q) = —^/C^^(q) for k > q and F k (q) = —l/(C k \q) — w(q)) for k < q which allows us 
to rewrite the second contribution to ()2.64|) as 



II = -k 



D-l 



dn ' k Cjpfa k, -p - fc)Cf (-p, -k,p + k)F k (k + p) - M ~ )M !;l - ■ (2.65) 



(2tt 



i-D 



l + T k {k)w{k) 



where we made use of lemma (I2.58J) . Collecting the intermediate results ()2.63j) and (j2.65|) 
we obtain finally 



dkC™ (p) 



= k - 1 

+ V- 1 

2 



^C { k 3 \p,k,-p-k) 2 'F k (k + p) 



(2< 



l + F k {k)w(k) 



dVLt 



C { k\p, ~P,k,-k) 



(27T 



T k {k) 2 w{k) 
1 + F k (k)w(k) 



[2.66) 
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in agreement with Parola and Reatto's result. Obviously, flow equations for the higher order 
Cfc"^ {n > 3) will also be ill-defined in the scaling regime k —>■ 0. This observation raises the 
problem of the validity of HRT in this regime and could be related with the singularities of 
high (~ infinite) compressibility states of HRT discussed recently by Reiner—. 

III. THE RG IN KSSHE FIELD THEORY 

These past few years two rigorous scalar field theories have been developed to describe 
simple fluids: the collective variable (CV) method introduced by the Ukrainian schoo l 47 ! 48 ' 49 
and the more recent KSSHE field theor y 17 ' 18 ' 19 . The CV theory involves two scalar fields 
with simple physical interpretations while only one field is involved in KSSHE theory with 
however a quite obscure physical signification. Both theory were shown to be equivalent 
order by order in the loop-expansion recentl y 20 ' 21 . In order to relate the RG equations of 
Section II with RG equations introduced these last past years in field theory and notably 
in Wetterich works 8 the KSSHE framework is at first sight more handy and will be adopted 
here. In this framework one can rewrites the GCPF (|2.2j) of a fluid as a functional integral 
thanks to a Hubbard-Stratonovic h 33 ' 34 transformation^ 9 -. For attractive pair potentials one 
has exactly (i.e. without unspecified normalization constant) 



where ip is a real random scalar field, w -1 is the inverse of w (in the sense of operators, i.e. 
w(l, 3) • w~ 1 (3, 2) = 5(1, 2)), and Shs \p + denotes the GCPF of bare hard spheres in the 
presence of a local chemical potential V(x) + tp(x). H can thus be written as a Gaussian 
average (. . .) w of covariance w and we have noted by M w the normalization constant 



The functional integrals which enter equations (|3.1|) and ()3.2|) can be given a precise meaning 
in the case where the domain Q is a cube of side L with periodic boundary conditions, in 
this case the measure Vtp reads as^ 




(H HS [v + ip\) 



w i 



(3.1) 




(3.2) 




dtpg 



(3.3a) 




2 ddtifq dQ(pq for q ^ . 



(3.3b) 
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where A = (2tt/L) IP (Z ring of integers) is the reciprocal cubic lattice. Note that the 
reality of <p implies that, for q ^ 0, (p- q = (p*, where the star means complex conjugation. 
With this slick normalization of the functional measure one has exactly 

M w = exp(^ J\nw(q)J ( 3 - 4 ) 

in the limit of large systems (L — > oo). Unless some UV cut-off k max is introduced the 
normalization constant M w is infinite, but even in the limit k max — ► oo the ratio which 
defines the GCPF in equation (|3.1|) will remain finite and independent of the cut-off for 
sufficiently large k max . 

It is instructive to rewrite equation (j3.1|) as 



E[u\ = U~ l j Vp exp (-W K [i/,^]) , 



^kK^] = ^(<f\w V) - ln -hs [v + (p] , (3.5) 

where Hk [v, <p] is the tree-level action of the KSSHE field theory since it follows from our 
dicussions of Section II that at scale-k of the RG group we can rewrite Ek [u] as 

Sfc M = j V<p exp (-H K W, <p\-\ {^RkWVj ■ (3-6) 

Indeed (1, 2) = it? -1 (l,2)+i?jfc(l,2) as can be inferred from the definition (|2.2U|) of the 
pair interaction Wk at scale-k. The presence of the normalization constant M Wk in the above 
equation is mandatory in order to obtain the exact flow of lnS^] and (3Ak [p], otherwise 
they are defined only up to an additional constant (i.e. which depends only on scale-k and 
not on the density p). 

Connoisseurs of NPRG theory will have noticed that the family of models introduced 
at Equation (|3.6J) is very similar to that considered by Wetterich^^ for generic scalar field 
theory. His idea was to build a one-parameter family of models for which a momentum- 
dependant mass term | j Rk(q)ip q <p_ q has been added to the original Hamiltonian (here 
Hk [v, <p}) in order to decouple the slow modes by giving them a large mass 8 . Thus equation 
()3.6|) could have been our starting point to obtain the RG flow of ln and f3Ak- By passing, 
the flow equation ()2.30j) for lnSj, can be obtained in the framework of KSSHE field theory 
without any recourse to the GC formalism. This somehow redundant material is left as an 
exercise for the reader. 
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A more detailed comparison of the NPRG of section II with the NPRG for generic (or 
canonical) field theory is worthwile but not so trivial. We have pointed out the similarities 
of the two RG constructions but there is a fundamental difference. The KSSHE theory is 
not a canonical scalar field theory in the sense that the coupling between the external source 
v(x) and the field <p(x) is non-linear as apparent in equation ([3.5)1 . In order to establish the 
exact link between the RG group of Section II and the NPRG of Wetterich we need first to 
build a canonical field theory for simple fluids, this is the aim of next section. 

IV. THE RG IN A CANONICAL FIELD THEORY OF SIMPLE FLUIDS 
A. A canonical field theory for simple fluids 

The starting point will be the equation ()3.5|) for 5 [u] that we rewrite with a slight change 
of notations (x instead of <p to avoid further confusions) as 

SM = A C 1 / V X exp^-i<x|^ 1 |x> + lnS H s^+x]) (4-1) 
Let us make now the change of variables 19 

V+x = Vo + V (4.2) 

where vq is an arbitrary constant, hence (p = \ + Av with Av = v + v$ — vq. With this trick 
one is led to a canonical field theory defined by its //-independent Hamiltonian 

W [<p\ = g - ln5 HS N + V?] , (4-3) 

and its partition function 

~* [J] = AC 1 J V<p exp {-U [cp] + (J\<p)) , (4.4) 
where the external source 

J(l) = ty- 1 (l,2) • Ap(2) . (4.5) 

Denoting by W* [J] = In 2* [J] and W [v] = In H [v] the grand potential (up to a multiplica- 
tive —ksT) we have of course 

W* [J] = W[u] + - (Az/|w _1 |Av) VJ, v given by Av = w ■ J , (4.6a) 

W [v\ = W* [J] - \ (J\w\J) Wu, J given by J = w' 1 ■ Av . (4.6b) 
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The above relations may be used to establish the relations between the density correlation 
functions T and the truncated correlation functions of the field (p, T (l, . . . , n) = 
5W* [J] /5J(1) . . . 5J{n). Noting that it follows from equation (@~5J) that 5 [. . .] /5J(1) = 
w(l, 2) ■ 5 [. . .] /Si/(2) on obtains readily from equations (J4.6|) that 

$(1) = (<p(l)> = w(l, 2) • p(2) + Ai/(1) , (4.7a) 
G£> T (1,2) = w(l,l')w(2,2')G^ T (l',2') + w(l,2) (4.7b) 
G$> T (l, ...,n) = w(l, l') . • • w(n, n) G {n) T (l', . . . , n) Vn > 3 . (4.7c) 

Since the KSSHE field \ and (p differ by a constant their truncated correlation functions 
coincide for n > 2, and, for n = 1 of course $(1) = (x(1))k + ^ u (^) fr° m which the relation 
p(l) = u; _1 (l, 2) (x(2)) K follows (the subscript [. . .] K means that the statistical average must 
be computed in the framework of KSSHE field theory). Therefore equations ()4.7|) establish 
the link between the GC formalism and the KSSHE field theory as well. 

We turn now our attention to the Gibbs free energy T [$] which will be defined as usual 
as the Legendre transform of the convex functional W* [Jp^ : 

(V0 G J 7 ) r [$] = sup {( J|$) - W* [J]} , 

Jej 

= (J*\§) - W*[J*] (4.8) 

where J* is the unique solution of the implicit equation 

5W* [J] 



= $(1) , (4.9) 

j=j* 



SJ(1 

and T and J are the convex sets of fields and sources respectively. One has then of course 
J(l) = 5T [$]/5$(l). 

We will now establish the link between the functionals $ i — > T [$] and p i — > (3 A [p\. To 
unclutter the notations we will consider rather 

PA\p]±/3A\p] + (^s-^o|p> , (4.10) 

so that 0A[p] = sup pen {(p\Au) - W[v)} and Ai/(1) = 50A [p] /5p(l) (since /3^4[p] and 
(3 A [p] differ by a linear functional, /L4. [p] is also a convex functional of the local density). 
In the one hand the relation (|4.7aJ) between the order parameters p and $ can be rewritten 
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as 

5A[p] 



$(1) = w(l,2)-p(2) + 



5p(l) 



p(l) = w-\l,2) -$(2) - S|, (4.11) 

and, on the other hand, replacing [J] by its expression ()4.6|) in term of the grand potential 
W [v] yields the useful identity 

(3A[p]=rm -i(J|Ai/>, (4.12) 
from which we deduce that either 

( v, eK )^w = r[*]-iraK/ r l*l 



2 \ 5$ 1 1 5$ 



d>sol.p(l) = W - 1 (1,2).$(2)-^M, (4.13) 
where "sol." means "solution of, or 

( v*^)r W = MM 

p sol. 0(1) = „(l,2).p(2) + ^M. (4.14) 

Finally the vertex functions of the canonical theory will be defined as usual^ 

T (n) m (1 n) = (4 15) 

It is also possible to establish the general relations between the [$] (1, . . . , n) and the 
direct correlation functions of the fluid C^ n > [p] (1, . . . , n) but they look rather awkward. By 
the way, only the cases n = 1, 2 will be useful, we quote them without the proofs which are 
elementary 

CW(1) = • U ;- 1 (l,2)-rW(2) , 
C (2) (l, 3) ■ T (2) (3, 2) = w-\l, 3) ■ C (2) (3, 2) + w(l, 3) ■ T (2) (3, 2) . (4.16) 



B. The RG flow in the canonical theory 

Following Berges et a/.— one defines 

[A = J Vip exp f-H [<p] - l - (y\R k \y) + (J\^ , (4.17) 
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where H, [cp] is the Hamiltonian ()4.3|) of the canonical theory of the fluid, Rk is the regulator 
defined in Section III B 11 and M Wk the normalization constant (|3.2|) (with the replacement 
w — > Wfc). The normalization constant M Wk is rarely considered in the literature on statistical 
field theory but here, it ensures that S£ [J] is exactly the partition function of the canonical 
field theory of a HS fluid with Wk pair interactions the bare action of which is given by 
H k [<p] = H[<p] + (<p\R k \<p) /2. 

The effective average action r k [<&] is defined as^ 

r' k m = sup {(jm-w* , (4.i8) 

(where W k [J] = lnH£) by the relation 

r;[$]=r fc [$] + i($| J R fc |$) . (4.19) 

Note that T' k [<&] is convex while [$] is not. 

In the limit k — > 0, — > and thus r^ =0 = Tk=o = T. In order to obtain the k — > oo 
limit of Tfc we start from its functional implicit expression 8 

exp (-r fe [$]) = AC, 1 J V X exp f« [$ + X ] + (xl^^) " \ (xi^lx)) • (4-20) 

For — > oo, — > oo and therefore ~ R k . The point is that we have exactly, in the 
limit Rk — ■> oo, without unspecified multiplicative constant, 

' e W (-Ux\R k \x) ) =*Ax\ , (4.21) 



where the functional Dirac mass 5jr [x] is given, with Wegner's normalizatio n 2 ! 20 ' 21 by 

Sjr ix] = V2W 5 (xo) II \ W 5 5 ^)\ > ( 4 - 22 ) 

as a direct calculation will show. Reporting equation ()4.21j) in ()4.20j) we obtain Tk=oo [$] = 

The RG flow of Tk [$] is given by the equation 8 

d k r k m = \J dld2d k R k (l,2)[R k + T k 2) y\l,2) 

- l -Jdld2 d k R k (l, 2)w k (l, 2) , (4.23) 
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(2) 

where T K k } is one of the effective vertex functions 

if > m (1 n) = 6Tk [$] (4 24) 

The first contribution to the right hand side of equation ()4.23|) seems to have been obtained 
for the first time by Nicoll^ and rediscovered later by Wetterictf while the second one is 
not considered in general and stems from the trivial flow of the normalization constant M Wk . 
Indeed it follows from equation (|3.4jl that 



d k \nN Wh = d k J\nw k {q) 

Rk(q)w k (q) 



dld2d k R k (l,2)w k (l,2) . (4.25) 



_V 
~ ~ ~2 
1 

~ ~2 

The equations for the flows d k T^ [<3>] (1, . . . ,n) of the running vertex are obtained by 
taking the functional derivatives of equation ()4.23j) with respect to $(x)&. They are 
obviously identical to those derived in section pi Cj) for the effective correlation func- 



tions under the replacements Cl — > — 1^ and d k w k — > d k R k for the vertex and 



J k 1 k 



G&> T = - (cf } + w k - u>} — > Gf' T = {if } + R k \ for the propagators. 

Until now we have only considered the first step of Wilson program^, i.e. the coarse 
graining operation. In order to get fixed points and apply the whole machinery of the RG 
group theory to the study of critical phenomena we have to consider now the second step, 
i.e. a rescaling of distances, wave vectors and fields. We refer the reader to the abundant 
literature on this poin t 1 ! 5 ' 8 ! 9 ' 25 , and only quote the RG equation obtained after rescaling 

d t T t = g^Y t + l - J dl d2 d t Rt{l,2) {Rt + r^y 1 (1,2) , (4.26) 

where t = ln(k/k max ) is the RG "time" and the expression of the dilatation operator may 
be written as 5 - 

Sdii = / q- dfi q _£- + (£>- D # ) / , (4.27) 

Jq 5§ q Jq 5$ q 

where D§ = (D — 2 + rj)/2 is the anomalous dimension of the field $ (recall that Fisher's 
exponent rj is a positive constant defined with respect to a non trivial fixed point of the 
flow ()4.26j) ). Note that we have discarded the trivial flow of the renormalization constant 
A/t from equation ()4.26|) to display an uncluttered formula. The flow of (3A k can be treated 
in a similar way 25 . 
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We want to point out that by contrast with the flow equations obtained in the GC or 
KSSHE theories, that obtained in the canonical theory of liquids is independent of the pair 
interaction w(r). All simple fluids are thus described by the same flow, they differ only by 
the initial conditions at k = oo which prescribe the microscopic model, i.e the pair potential 
and the reference system (something else than hard spheres can be chosen). This property 
together with the existence of fixed points and a group law for the flow (|4.26|) yield the 
universality of critical properties of all simple fluids with short range interactions. Their 
universality class is of course that of the Ising model. 

C. The relation between the flows of F k [<£] and f3Ak [p] 

It transpires from the previous sections that the flows of Tk [$] and PAk [p] must be 
related in some way since both describe the same physics. We first note that the k-system 
being a HS fluid with pair interactions Wk(r) the relations established in Section (jlV Aj) 
between the "true" action T' k and the "true" GC free energy [3A k can be applied at scale-k 
and one has for instance (cf equation ()4.13|1 ) 

PA k [ P ] = r fc [$]--(— - 



2 \ 5$ 11 5$ 

$sol.p(l) = ^(1,2) -$(2) -^21, (4.28) 

where 

PA k [p]=PA k [p] + (u s , k -u \p) ■ (4.29) 

Recall that PA' k [p] and PA k [p] are related through equation (j2.24jl . Now we define 

P~A k [p] 4 pA k [p] + (ua - u \p) (4.30) 

so that equation (|2.24|) can be rewritten as 

PA k [p] = PA k [p] + ~ (p\w k - w\p) . (4.31) 

Reexpressing equations (|4.28jl in terms of Tk and Ak thanks to equations (j4.19|) and (j4.31|) 
yields 

PA k [ P ] = r k m--(^-±\w\ [ 



2\ 5$ 11 5$ 
$ sol. p(l) = w-\l, 2) • $(2) - . (4.32) 
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This remarkable result shows that the functional relations between the running Tk and (3Ak 
as well as those between the order parameters p and $ are in some sense independent of 
scale-k. Proceeding in the same way with equation (|4.14j) one finds the dual relation 

r t [*]=/3A M + i/»lk-' l ^W 



2\ dp 1 1 Sp 
psol.$(l) = Hl,2).p(2) + ^M. (4.33) 

Note that the mapping p > $ is a one to one correspondence since we have, for instance from 
equation (jOSjl . that 5$(l)/5p(2) = w(l, 2) -Cf 2), which is a positive definite operator 
since it follows from equation ()2.28|) that w(q) — C^\q) > Wk{q) > 0. For homogeneous 
systems the order parameter <3> is therefore an increasing function of the density p. As a 
consequence of equations ()4.32|) and ()4.33|) we also have 

^M = W (12).^M (4 34) 

or equivalently (1) = — w(l,2) • 1^ (2). This relation between n = 1 vertex functions 
and direct correlation functions happens once again to be an identity independent of scale-k. 

We are now in position to establish the relation between the two flows dk(3Ak [p] and 
dkTk [$]. Taking the partial derivative of equation (|4.32jl with respect to k at fixed p, one 
finds that 

/5f3A k i _ii „ S0Ak 



dk(3A k [p} = dkTkin.-^j^^dk- 



5 P _ 

= 9 fc r fc [$] + ^ fe — | W . — (4.35) 
where the last line follows from 

9k W) = (436) 
as can be inferred from equation ()4.33|) . Note that the second term in the right hand side 
of the last line of equation ()4.35j) vanishes by virtue of equation (|4.34j) . Since dk/3Ak [p] = 
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dk/3Ak [p] (cf equation (|4.H()jl ) we have finally 

d k (3A k [p] = d k T k [$] , 
(Vp), $sol. P (l) = ^(1,2) ■ $(2) -^M 



or 

(V$), p sol. $(1) = u>(l,2) • p(2) + + i/ s - u • (4.37) 

Therefore the two flows coincide, moreover it is not difficult to check that the initial condi- 
tions PAoo [p] = PAuf [p] and Too [$] = T~t W\ are of course also related by relations ()4.37j) 
as a short calculation will show. 



D. Additional Comments 

We complete the results of Section (jlV CJ) with a study of correlation functions at scale-k. 
The relations ()4.7j) are of course valid for the k-systems, so that we have 

G®?(1,2) = Wfe (l,l'K(2,2')Gi 2)T (l',2') + w k (l,2) (4.38a) 
G { ;\ T (l,...,n) = w k {l,l')...w k (n,n')G k n)T (l',...,n) Vn > 3 , (4.38b) 

where the G^ n l T are the correlation functions of the KSSHE field at scale-k. The relations 
(I4.38|) may be used to derive the identity of the flows (|4.37|) starting from equation (|2.32|) 
for d k (3A k by replacing G k 2 ^' T by G^' k T by its expression (|4.38aj) : one then obtains the 
expression ()4.23|) for d k T k (including the flow of the normalization M Wk ) as the reader will 
check easily. 

Concerning the vertex functions we have already given the relation between C k and 
rj c n_1 ' ) (cf equation (|4.34Jl ). For n = 2 and at scale-k we have from equation (J4.16j) 

q (2) (l, 3) ■ rl (2) (3, 2) = ^(l, 3) ■ C; (2) (3, 2) + ^(1, 3) • r'« (3, 2) , (4.39) 

from which it is easy infer that 

Cf(l, 3) • if >(3, 2) = ^(1, 3) • Cf (3, 2) + w(l, 3) • lf>(3, 2) . (4.40) 

It is likely -but we did not try to prove it- that the relations between the Cj^ and rl n ^ are 
also independent of the scale-k at any order 'n', i.e. they involve w rather than w k . 
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From this lengthy discussion on the identity of the flows dkfiAk and we can draw 
general conclusions which enlighten some obscure points related to the singular behavior of 
the HRT theory in the subcritical region. Known properties of the flow of I\ can be trans- 
posed to the flow of j3 Ak thanks to the exact equations (|4.32j) and (|4.37j) . For instance, it 
is known from theoretical and numerical studies of many different models treated in various 
approximation schemes that, in general, d 2 Tk/d<& 2 is discontinuous at the coexistence curve 
below T c , jumping from a finite value in the ordered phase to zero in the inner region of the 
coexistence curve. In other words, the susceptibility jumps from a finite value to infinity. 
This property survives the one-to-one mapping ()4.32j) yielding the conclusion that, for simple 
fluids, the second derivative of the free energy d 2 fk/dp 2 will also exhibits a finite disconti- 
nuity on the binodal, being zero in the two phase region, as a short calculation will show. 
By contrast the HRT yields a continuous d 2 fk/dp 2 25 i 26 i 27 i 28 and thus, in this framework, 
the binodal and the spinodal coincide. It is likely that the ultra sharp cut-off limit needed 
to extract HRT from our smooth cut-off equations enforces the continuity of d 2 fk/dp 2 . One 
could argue that this flaw of the theory could also be due to the approximations injected 
into the exact HRT flow hierarchy in order to solve it. We do not think so and support the 
former hypothesis basing our point of view i) on the fact that the discontinuity of d 2 Ti i /d^ 2 
at subcritical temperatures is observed whatever the approximation retained^ ii) on the anal- 
ysis of Section pi Cj) showing that the HRT equations for the two-body and higher-order 
direct correlation functions are ill-defined in the scaling regime k — > 0. If this is true then 
there is no way to improve HRT even by incorporating sophisticated closure relations such 
as those discussed in reference^. Unfortunately, in his own smooth cut-off formulation of 
HRT (which however cannot reduce to our formulation in any way, as discussed in Section 
(II)) Parola does not discuss this point^. We cannot help raising the doubt: is HRT really 
intrinsically correct? 



V. RG ANALYSIS OF KAC MODEL 

As an illustration of the issues of previous sections we study here the RG flow of Kac 
model which can done analytically. 
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A. Generalities 



In the Kac model the hard spheres interact via a long range pair potential w(r) obtained 
as the limi t 22 ' 31 ' 37 ' 38 ' 39 

w(r) = lim a D &(ar) at fixed volume V , (5-1) 

where $(r) is a short-range function, distances are measured in unit of a and the inverse 
range of the interaction a > is dimensionless. We assume <3>(fc) > and $(0) finite and 
note that w(k) = §(k/a). It is important to emphasize that the limit a — > must be 
taken before taking the thermodynamic limit, i.e. at fixed finite volume V. Assuming PBC 
conditions the potential energy E P (C) of any GC configuration C will be 

-E p {C)/k B T = w(k)p(k)K-k) ! (5.2) 

fceA 

where p(k) is the Fourier transform of the microscopic density 1)2.1)1 . For a given volume V, 
w(k) vanishes for any k ^ of the Fourier space A in the limit a — > while w(0) = $(0). 
Therefore the potential energy E P {C) of any GC configuration C will be 

-E p (C)/k B T = ^p(0) 2 $(0) 

= ^P 2 , (5-3) 

where we have changed our notations and defined the reduce inverse temperature as f3 = $(0) 
and p = N/V is the number density in configuration C. The GCPF of the model is then 
given (for a homogeneous system) by 

S£ ac [v] ~ J™ dp exp (-V (-up + f m (p) ~ \PP 2 ^J \ , (5.4) 

where we have replaced the sum over in equation ()2.2)) by an integral over the density 
which is valid for large systems with no consequences on the thermodynamic limit and /hs(p) 
is the specific Helmoltz free energy of the HS fluid at density p. Note also that, in the limit 
a — > 0, one has V — v since the self-energy w(0) vanishes. Only now we take the limit 
V — > oo with the result 



lim V~ L \nEy[i/} = min{£(i/,p)} , (5.5) 
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where 

£(v>p) ~ /mf(p,/?) - pv , 
/mf(p,/3) ^ /hs(p)-§ P 2 - (5.6) 

The function £(z/, p) of equation ()5.5|) plays therefore the role of a Landau function^. Note 
that /mf(Pj/5) is precisely the MF free energy as defined at equation (j2.14j) . The HS free 
energy /hs(p) which enters the expression of C(u, p) is known analytically only in dimensions 
D = 1 and D = oo, however in the case D = 3 -to which we stick to from now- good 
approximations are known, for instance the Carnahan-Starling (CS) approximation^ for 
which we have 

/ H s(p)=p(lnpA 3 -l) + 6 ^ 4 ~y , (5.7) 

where A is de Broglie wavelength and rj = 7rper 3 /6 is the packing fraction. The CS approxi- 
mation describes a single fluid phase from rj = up to the unphysical r\ — 1 and the function 
of equation ()5.7|) is strictly convex, i.e. its second derivative / HS (p) > for all < rj < 1. 
In fact the function p i— > f^(p) > is itself convex with a unique positive minimum at 
p a 3 = 0.2491295 . . . with (3 C = /4s (Po) = 11.101579 . . . i&. It follows from these remarks 
that for f3 < f3 c the Landau function p i— > £(v, p) is convex and coercive for any v (since 
<9 2 /mf(p, P) /dp 2 > (P c — 0) for all p) and that the minimum of C(u, p) in equation (j5.5j) 
exists and is unique. Conversely for (3 > j3 c , C(v, p) admits in general two minima and there 
is only but one value of the chemical potential v = z/ coex (/5) for which the two minima p g ((3) 
and pi((3) (pi > p g ) which describes the gas and the liquid phases respectively are such that 
^(^cocx,P s ) = £(^cocx,pO = -flPcaexifl) where /3P coe x is the pressure at coexistence. At (3 C 
the two densities p g and pi merge to po and thus f3 c is the inverse critical temperature of 
Kac model and po its critical density. Note that for f3 > (3 C the Kac pressure has a kink at 
v = z^cocx with two distinct slopes p g and p\. To this kink corresponds a segment for the 
Legendre transform, i.e. the GC free energy^. One thus have, for j3 < f3 c 

/ Kac (p,/3) = /mf(p,/5) , (5.8) 

and for ft > f3 c 

/Kac(p,/3) = /mf(p,/?) for p> Pi(ft) or p < p g (j3) , 

= fuF^ft)-^ 1 + /mf(pz,/?)^^ for Pg < p < pi , (5.9) 
Pg ~ Pi Pi- Pg 
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FIG. 1: Flow of f Y {p) - v mex {^)p at /3 = 13 > f3 c . At Y = f3 (thick solid line) f Y = /mf, at Y = 
(thick dashed line) /y = /kac- As Y decreases (dashed lines from top to bottom) fy becomes 



convex. 



where p g (P), Pi{0), and z/ coex (/3) are the unique solution of the system 

fMF{Pg,P) — V coex p g = fuFiPliP) — ^cooxPi , 

dfMF(Pg,(3)/dp g = df M F(pi,P)/dpi = z/ cocx . 



(5.10) 



B. NPRG for Kac model 



Our starting point for the NPRG analysis of Kac model will be equation ()2.37|) that we 
rewrite as 



df k (p)/dk 2 



dw k (q)/dk' 2 



(5.11) 



2j qC P(q)+w k (q)-w(q) ' 
We do not aim at mathematical rigor and will solve the flow by hand-waving arguments. 
Taking the Kac limit a — > at fixed V in the RG framework means that at any scale k ^ 
one has k 3> a and g>«, therefore 



w(q)~(l-e a (q))P, 



(5.12) 
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where 6 a (q) is a smoothened step function equal to 1 for q > a and to for q < a. Whatever 
the form of the cut-off function Rk(q) one thus has 

w k ( q ) ~ (i - e a (g)) 1 _f k2j3 , 

d k 2 Wk (q) ~ (l-9 a (g)) 2 , (5.13) 

(1 + fc 2 /?) 

since Rk{q) ~ & 2 for g — > (cf the conditions (|2.21j) . Because of the presence of the 
(1 — G Q (q)) in the numerator of the right hand side of equation ()5.11j) one can replace the 
function C^\q) at the denominator by its value at q = 0, i.e. —d 2 fk(p)/dp 2 . Therefore for 
a ^ k one has the simplified flow equation 

4vr , (3 2 



df k (p)/dk 2 = —a* j-Z , (5.14) 

6 (1 + k*fi) (d 2 f k /dp 2 + gfa 

which can be rewritten, with the change of variable 

k 2 f3 2 
l + k 2 (3 

as 



Y = 7^775, (5-15) 



df Y (p)/dY [d 2 f Y (p)/dp 2 + Y]= —a 3 , (5.16) 

where fy(p) = fk{p)- One can take safely the limit a — > of the equation above which 
yields 

df Y (p)/dY [d 2 f Y (p)/dp 2 + Y] = . (5.17) 

This somehow strange partial differential equation must be solved with the initial conditions 
(|2.29J1 . i.e at k = oo 4=> Y = [3 one has f Y=/3 (p) = /mf(p, P) = /hs(p) — Pp 2 /2- us define 
X = Y + d 2 fY(p)/dp 2 . We first remark that, at the beginning of the flow, i.e. for Y = /3, 
the initial conditions are such that that X = /jj s (p) > (3 C = /hs (Po) > for all p. It follows 
then from the flow equation (|5.17|) that dfy (p) / dY = and therefore the flow does not start 

(2) 

immediately. In this initial stage of immobility X = [p) — (3 + Y > (/3 C — /3 + K ) for all 
p, while df Y (p)/dY = if /5 C - /3 + F > 0. 

Now, if /3 < /3 C then X remains strictly positive for all p when Y varies from Y = /3 
(k = oo) to Y = (A; = 0). Therefore dfy{p)/dY remains zero and the flow never starts. 
We thus have fy=o{p) = /mf(p,/3) = /Kac(p,/3) as expected. 

The case /3 > /3 C is more interesting. For K = (3 — (3 C the minimum of X(p) becomes 
equal to zero and the flow can start. If it actually does then it follows from ()5.17j) that 
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FIG. 2: Flow of dfY(p)/dp at (5 = 13 > (3 C . As F decreases (dashed lines from top to bottom) a 
plateau appears. 

Y + d 2 /y(p) / 'dp 2 = for < Y < Y in some domain of p. It seems reasonable to assume 
that, at given Y, this domain is some interval [p\{Y), p 2 {Y)}. Therefore 

f Y {p) = - l -Yp 2 + a{Y)p + b{Y) (5.18) 

for pi < p < p 2 . For p > p 2 or p < pi the flow has not yet started and fy remains equal to 
its initial value /mf(p, /?)■ In order to determine the unknown functions a{Y), b(Y), pi(Y) 
and p 2 {Y) we invoke the continuity of the functions p t— > /y(p) and p \— > dfy /(p)/dp which 
can be inferred from the general flow equation ( I5.11|) . Hence we have VF G [0, Y] 

fMPuP) = ~Yp 2 + a(Y) Pi + b(Y) (i = l,2), 
df M F(Pi,/3)/d Pl = -Y Pi + a(Y) (i = 1,2). (5.19) 

In order to solve equations (|5.19j) we remark that the continuity conditions 1)5.19)1 can be 
reformulated as /mf(p«, A) = a(Y)pi + b(Y) and dfupiPi, A) /dp = a(Y) with A = (3 — Y. 
These relations are precisely those which determine the coexistence curve of the Kac model 
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at the inverse temperature A (cf equations ()5.10|1 ). Therefore 



Pl (Y) = p g (A) , 

p 2 (Y) = Pl (A) . (5.20) 

In other words, the flows of pi(Y) and p2(Y) are given by the two branches of the binodal 
of the Kac model at inverse temperature A = (3 — Y. The starting point of the flow Y = Y 
corresponds to A = /3 C , i.e to the inverse critical temperature of Kac model. Moreover 
a(Y) = z/ coex (A) coincides with the coexistence chemical potential of the Kac model at 
inverse temperature A and b(Y) = —/3P coex (A) with minus the pressure at coexistence. Our 
final result can be summarized as 

fy(p) = /mf(a/3) for p < p g (A) or p > p t (A) , 

= --Yp 2 + iw(A)p - /3P cocx (A) for Pg (A) <p< pi{A) . (5.21) 

At the end of the flow, i.e. for Y — 0, one has A = f3 and thus fy(p) = /KacG°> P) as given 
by equations ()5.8|) and (|5.9p . 

To conclude this section let us have a look on the flow of the "true" free energy f' k (p). It 
follows from equation (J2.24j) that 

fl(p) = f k (p) + lp 2 (w(0)-w k (0)) , 

= f k (p) + l -Yp\ (5.22) 

from which we deduce that f Y (p) = f'k(p) = /Kac(p, A) which is an obvious result since 
the k-system is nothing but a Kac model at inverse temperature A as can be inferred from 
equations (j5.12j) and (J5.13)) . This could have been our starting point to solve the flow. This 
allows to find the exact solution of the flow without solving explicitely equation (|5.11j) . but 
it is less funny. 

To illustrate our result we have displayed the flow of the running free energy and its 
two first derivatives with respect to p. Figure 1 shows the flow of /y(p) — t / CO ex(/3)p at a 
temperature T* = 1/(3 below T c as Y varies from Y = (3 to Y = 0. This demonstrates 
the approach of convexity in the inner region when Y —>■ 0. Figure 2 displays the flow 
of dfy(p)/dp at the same subcritical temperature, demonstrating the slow construction of 
the plateau. Finally, in Figure 3 we displays the flow of d 2 fy(p)/dp 2 which exhibits a 
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FIG. 3: Flow of 9 2 /y(p)/9p 2 at [3 = 13 > /3 C . The compressibility is discontinuous as soon as the 
flow starts. 

finite discontinuity at the points pi(Y) = p g (A) and P2OO = Pi(A). Since at Y — 0, the 
isothermal compressibility xt given by /^x^ 1 = p~ 2 d 2 f Y =o(p)/dp 2 , xt is finite on the binodal 
although infinite in the coexistence region. 

We end this section by pointing out that the HRT flow cannot be solved conveniently for 
the Kac model, since at any scale fc^Owe have Wk(q) = Wq, as a consequence of the ultra 
sharp cut-off. In this case the flow starts and stops exactly at k — with an uncontrollable 
singularity. 

VI. CONCLUSION 

In this paper we have developed a smooth cut-off NPRG theory for simple fluids. Two 
equivalent versions of the theory were formulated, the first one in the GC ensemble, with 
an immediate transcription in the framework of KSSHE field theory, the second one in the 
framework of a generic scalar field theory (with no Z2 symmetry). The exact mapping 
between these two versions was derived in detail. The HRT is obtained as a highly singular 
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limit of regular smooth cut-off NPRG equations. The inability of HRT to distinguish between 
the spinodal and binodal curves is probably a consequence of this limit which is only partly 
controllable (the flow equations for the with n > 2 are ill-defined in the scaling limit 
k — *■ 0). This shortcoming is absent of the present theory. 

Our developments where highly abstract and formal and we did not say anything on 
the strategy to adopt to solve the NPRG flow equations so some comments on this issue 
are probably wellcome. There are mainly two ways to solve the NPRG equations. A first 
method is to restrict the functional form of the effective average action Tfc[$] (or /3.4.fc[p]) 
on a restricted functional sub-space. For instance the simplest ansatz is the so-called local 
potential (LPA) approximation^^ 



where £/& (($(x))) is a function of the classical field. The LPA ansatz allready gives non 
trivial critical exponents (except i] = 0) and is quite easy to solve numerically^. Allready 
the improvement 1/2 — > in equation ()6.1|) which includes the renormalization of the 
wave function and yields rj ^ makes the results more satisfactory. It must be precised 
first that most of the numerical studies of LPA or improved LPA equations were applied to 
systems with Z2 symmetry, i.e. with no direct application to liquids, see however reference^ 
where symmetric and antisymmetric correction-to-scaling exponents are studied in detail. 
Secondly, these ansatz (including those with more terms in the gradient expansion) allow 
to estimate the low-q behavior of all the vertex but do not give their correct behavior at 
short distances in direct space. Such methods will not fit with the dogma of the theory 
of liquids according to which g{r) should be zero in the core. In liquids the equation of 
state (i.e. the virial pressure) results from a delicate balance between energetic and entropic 
contributions (i.e. the value of the g(r) at contact, that is at r = cr). If applied to liquid 
theory, approximations such as LPA or gradient expansions although giving good estimations 
of critical exponents will miss this point and would yield a poor estimation of non universal 
critical properties. This is why Reatto and coworkers considered another strategy to solve 
their HRT hierarchy. It consists in closing the tower of equations for the . Contrary to 
the first strategy high order are ignored but the structure of C^\r) can be described 
with good accuracy for all "r". In particular one can take advantage of the arbitrariness of 
Wfc(r) in the core to enforce gfc(r) = for < r < a at each scale-k. The most recent HRT 




(6.1) 
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studie a 27 i 28 i 29 i 51 consider a closure in the spirit of the Optimized RPA (ORPA) 52 

C?\r) = C£l(r) +I k w(r) + K k (r), 
K k {r) = J2 -r> n (r) , (6.2) 

n=l 

where the u n (r) are taken to provide a basis for a suitable function space over [0, a}. Coef- 
ficient 7° is chosen to enforce thermodynamic consistency (i.e (0) = —d 2 fk/dp 2 ) while 
the n max and the 7^ will be determined in such a way that g(r) is zero in the core with 
sufficient accuracy. It seems that this numerical scheme could be adapted to our NPRG 
flow equations, but other schemes proposed recently for solving Wetterich flow should also 
be considered^^i. We hope to report progresses on these points in future works. 
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